Single-fly assemblies fill major phylogenomic gaps across the Drosophilidae Tree of Life.

Long-read sequencing is driving rapid progress in genome assembly across all major groups of life, including species of the family Drosophilidae, a longtime model system for genetics, genomics, and evolution. We previously developed a cost-effective hybrid Oxford Nanopore (ONT) long-read and Illumina short-read sequencing approach and used it to assemble 101 drosophilid genomes from laboratory cultures, greatly increasing the number of genome assemblies for this taxonomic group. The next major challenge is to address the laboratory culture bias in taxon sampling by sequencing genomes of species that cannot easily be reared in the lab. Here, we build upon our previous methods to perform amplification-free ONT sequencing of single wild flies obtained either directly from the field or from ethanol-preserved specimens in museum collections, greatly improving the representation of lesser studied drosophilid taxa in whole-genome data. Using Illumina Novaseq X Plus and ONT P2 sequencers with R10.4.1 chemistry, we set a new benchmark for inexpensive hybrid genome assembly at US $150 per genome while assembling genomes from as little as 35 ng of genomic DNA from a single fly. We present 183 new genome assemblies for 179 species as a resource for drosophilid systematics, phylogenetics, and comparative genomics. Of these genomes, 62 are from pooled lab strains and 121 from single adult flies. Despite the sample limitations of working with small insects, most single-fly diploid assemblies are comparable in contiguity (>1Mb contig N50), completeness (>98% complete dipteran BUSCOs), and accuracy (>QV40 genome-wide with ONT R10.4.1) to assemblies from inbred lines. We present a well-resolved multi-locus phylogeny for 360 drosophilid and 4 outgroup species encompassing all publicly available (as of August 2023) genomes for this group. Finally, we present a Progressive Cactus whole-genome, reference-free alignment built from a subset of 298 suitably high-quality drosophilid genomes. The new assemblies and alignment, along with updated laboratory protocols and computational pipelines, are released as an open resource and as a tool for studying evolution at the scale of an entire insect family.

The Transmission Patterns of the Endosymbiont Wolbachia within the Hawaiian Drosophilidae Adaptive Radiation.

The evolution of endosymbionts and their hosts can lead to highly dynamic interactions with varying fitness effects for both the endosymbiont and host species. Wolbachia, a ubiquitous endosymbiont of arthropods and nematodes, can have both beneficial and detrimental effects on host fitness. We documented the occurrence and patterns of transmission of Wolbachia within the Hawaiian Drosophilidae and examined the potential contributions of Wolbachia to the rapid diversification of their hosts. Screens for Wolbachia infections across a minimum of 140 species of Hawaiian Drosophila and Scaptomyza revealed species-level infections of 20.0%, and across all 399 samples, a general infection rate of 10.3%. Among the 44 Wolbachia strains we identified using a modified Wolbachia multi-locus strain typing scheme, 30 (68.18%) belonged to supergroup B, five (11.36%) belonged to supergroup A, and nine (20.45%) had alleles with conflicting supergroup assignments. Co-phylogenetic reconciliation analysis indicated that Wolbachia strain diversity within their endemic Hawaiian Drosophilidae hosts can be explained by vertical (e.g., co-speciation) and horizontal (e.g., host switch) modes of transmission. Results from stochastic character trait mapping suggest that horizontal transmission is associated with the preferred oviposition substrate of the host, but not the host's plant family or island of occurrence. For Hawaiian Drosophilid species of conservation concern, with 13 species listed as endangered and 1 listed as threatened, knowledge of Wolbachia strain types, infection status, and potential for superinfection could assist with conservation breeding programs designed to bolster population sizes, especially when wild populations are supplemented with laboratory-reared, translocated individuals. Future research aimed at improving the understanding of the mechanisms of Wolbachia transmission in nature, their impact on the host, and their role in host species formation may shed light on the influence of Wolbachia as an evolutionary driver, especially in Hawaiian ecosystems.

Bacterial and fungal components of the gut microbiome have distinct, sex-specific roles in Hawaiian Drosophila reproduction.

Gut microbiomes provide numerous physiological benefits for host animals. The role of bacterial members of microbiomes in host physiology is well-documented. However, much less is known about the contributions and interactions of fungal members of the microbiome even though fungi are significant components of many microbiomes, including those of humans and insects. Here, we used antibacterial and antifungal drugs to manipulate the gut microbiome of a Hawaiian picture-wing Drosophila species, D. grimshawi, and identified distinct, sex-specific roles for the bacteria and fungi in microbiome community stability and reproduction. Female oogenesis, fecundity and mating drive were significantly diminished when fungal communities were suppressed. By contrast, male fecundity was more strongly affected by bacterial but not fungal populations. For males and females, suppression of both bacteria and fungi severely reduced fecundity and altered fatty acid levels and composition, implicating the importance of interkingdom interactions on reproduction and lipid metabolism. Overall, our results reveal that bacteria and fungi have distinct, sexually-dimorphic effects on host physiology and interkingdom dynamics in the gut help to maintain microbiome community stability and enhance reproduction.

Phenotypic disruption of cuticular hydrocarbon production in hybrids between sympatric species of Hawaiian picture-wing Drosophila.

Interspecies hybrids can express phenotypic traits far outside the range of parental species. The atypical traits of hybrids provide insight into differences in the factors that regulate the expression of these traits in the parental species. In some cases, the unusual phenotypic traits of hybrids can lead to phenotypic dysfunction with hybrids experiencing reduced survival or reproduction. Cuticular hydrocarbons (CHCs) in insects are important phenotypic traits that serve several functions, including desiccation resistance and pheromones for mating. We used gas chromatography mass spectrometry to investigate the differences in CHC production between two closely related sympatric Hawaiian picture-wing Drosophila species, Drosophila heteroneura and D. silvestris, and their F1 and backcross hybrid offspring. CHC profiles differed between males of the two species, with substantial sexual dimorphism in D. silvestris but limited sexual dimorphism in D. heteroneura. Surprisingly, F1 hybrids did not produce three CHCs, and the abundances of several other CHCs occurred outside the ranges present in the two parental species. Backcross hybrids produced all CHCs with greater variation than observed in F1 or parental species. Overall, these results suggest that the production of CHCs was disrupted in F1 and backcross hybrids, which may have important consequences for their survival or reproduction.

Widespread introgression across a phylogeny of 155 Drosophila genomes.

Genome-scale sequence data have invigorated the study of hybridization and introgression, particularly in animals. However, outside of a few notable cases, we lack systematic tests for introgression at a larger phylogenetic scale across entire clades. Here, we leverage 155 genome assemblies from 149 species to generate a fossil-calibrated phylogeny and conduct multilocus tests for introgression across 9 monophyletic radiations within the genus Drosophila. Using complementary phylogenomic approaches, we identify widespread introgression across the evolutionary history of Drosophila. Mapping gene-tree discordance onto the phylogeny revealed that both ancient and recent introgression has occurred across most of the 9 clades that we examined. Our results provide the first evidence of introgression occurring across the evolutionary history of Drosophila and highlight the need to continue to study the evolutionary consequences of hybridization and introgression in this genus and across the tree of life.

Highly contiguous assemblies of 101 drosophilid genomes.

Over 100 years of studies in Drosophila melanogaster and related species in the genus Drosophila have facilitated key discoveries in genetics, genomics, and evolution. While high-quality genome assemblies exist for several species in this group, they only encompass a small fraction of the genus. Recent advances in long-read sequencing allow high-quality genome assemblies for tens or even hundreds of species to be efficiently generated. Here, we utilize Oxford Nanopore sequencing to build an open community resource of genome assemblies for 101 lines of 93 drosophilid species encompassing 14 species groups and 35 sub-groups. The genomes are highly contiguous and complete, with an average contig N50 of 10.5 Mb and greater than 97% BUSCO completeness in 97/101 assemblies. We show that Nanopore-based assemblies are highly accurate in coding regions, particularly with respect to coding insertions and deletions. These assemblies, along with a detailed laboratory protocol and assembly pipelines, are released as a public resource and will serve as a starting point for addressing broad questions of genetics, ecology, and evolution at the scale of hundreds of species.

Analysis of Genomic Sequence Data Reveals the Origin and Evolutionary Separation of Hawaiian Hoary Bat Populations.

We examine the genetic history and population status of Hawaiian hoary bats (Lasiurus semotus), the most isolated bats on Earth, and their relationship to northern hoary bats (Lasiurus cinereus), through whole-genome analysis of single-nucleotide polymorphisms mapped to a de novo-assembled reference genome. Profiles of genomic diversity and divergence indicate that Hawaiian hoary bats are distinct from northern hoary bats, and form a monophyletic group, indicating a single ancestral colonization event 1.34 Ma, followed by substantial divergence between islands beginning 0.51 Ma. Phylogenetic analysis indicates Maui is central to the radiation across the archipelago, with the southward expansion to Hawai'i and westward to O'ahu and Kaua'i. Because this endangered species is of conservation concern, a clearer understanding of the population genetic structure of this bat in the Hawaiian Islands is of timely importance.

A Whole-Genome Scan for Association with Invasion Success in the Fruit Fly Drosophila suzukii Using Contrasts of Allele Frequencies Corrected for Population Structure.

Evidence is accumulating that evolutionary changes are not only common during biological invasions but may also contribute directly to invasion success. The genomic basis of such changes is still largely unexplored. Yet, understanding the genomic response to invasion may help to predict the conditions under which invasiveness can be enhanced or suppressed. Here, we characterized the genome response of the spotted wing drosophila Drosophila suzukii during the worldwide invasion of this pest insect species, by conducting a genome-wide association study to identify genes involved in adaptive processes during invasion. Genomic data from 22 population samples were analyzed to detect genetic variants associated with the status (invasive versus native) of the sampled populations based on a newly developed statistic, we called C2, that contrasts allele frequencies corrected for population structure. We evaluated this new statistical framework using simulated data sets and implemented it in an upgraded version of the program BayPass. We identified a relatively small set of single-nucleotide polymorphisms that show a highly significant association with the invasive status of D. suzukii populations. In particular, two genes, RhoGEF64C and cpo, contained single-nucleotide polymorphisms significantly associated with the invasive status in the two separate main invasion routes of D. suzukii. Our methodological approaches can be applied to any other invasive species, and more generally to any evolutionary model for species characterized by nonequilibrium demographic conditions for which binary covariables of interest can be defined at the population level.

Varieties of the highly dispersible and hypervariable tree, Metrosideros polymorpha, differ in response to mechanical stress and light across a sharp ecotone.

PREMISE: The drivers of isolation between sympatric populations of long-lived and highly dispersible conspecific plants are not well understood. In the Hawaiian Islands, the landscape-dominant tree, Metrosideros polymorpha, displays extraordinary phenotypic differences among sympatric varieties despite high dispersibility of its pollen and seeds, thereby presenting a unique opportunity to investigate how disruptive selection alone can maintain incipient forms. Stenophyllous M. polymorpha var. newellii is a recently evolved tree endemic to the waterways of eastern Hawai'i Island that shows striking neutral genetic differentiation from its ancestor, wet-forest M. polymorpha var. glaberrima, despite sympatry of these forms. We looked for evidence for, and drivers of, differential local adaptation of these varieties across the range of M. polymorpha var. newellii. METHODS: For paired populations of these varieties, we compared seedling performance under contrasting light conditions and a strong water current characteristic of the riparian zone. We also conducted a reciprocal transplant experiment and contrasted adult leaf anatomy. RESULTS: Results suggest that the riparian zone is harsh and that selection involving the mechanical stress of rushing water, and secondarily, light, led to significant reciprocal immigrant inviability in adjacent forest and riparian environments. The strongest adaptive divergence between varieties was seen in leaves and seedlings from the site with the sharpest ecotone, coincident with the strongest genetic isolation of M. polymorpha var. newellii observed previously. CONCLUSIONS: These findings suggest that disruptive selection across a sharp ecotone contributes to the maintenance of an incipient riparian ecotype from within a continuous population of a long-lived and highly dispersible tree species.

Horizontal Transfer of Bacterial Cytolethal Distending Toxin B Genes to Insects.

Horizontal gene transfer events have played a major role in the evolution of microbial species, but their importance in animals is less clear. Here, we report horizontal gene transfer of cytolethal distending toxin B (cdtB), prokaryotic genes encoding eukaryote-targeting DNase I toxins, into the genomes of vinegar flies (Diptera: Drosophilidae) and aphids (Hemiptera: Aphididae). We found insect-encoded cdtB genes are most closely related to orthologs from bacteriophage that infect Candidatus Hamiltonella defensa, a bacterial mutualistic symbiont of aphids that confers resistance to parasitoid wasps. In drosophilids, cdtB orthologs are highly expressed during the parasitoid-prone larval stage and encode a protein with ancestral DNase activity. We show that cdtB has been domesticated by diverse insects and hypothesize that it functions in defense against their natural enemies.

Reproductive Capacity Evolves in Response to Ecology through Common Changes in Cell Number in Hawaiian Drosophila.

Lifetime reproductive capacity is a critical fitness component. In insects, female reproductive capacity is largely determined by the number of ovarioles, the egg-producing subunits of the ovary [e.g., 1]. Recent work has provided insights into ovariole number regulation in Drosophila melanogaster. However, whether mechanisms discovered under laboratory conditions explain evolutionary variation in natural populations is an outstanding question. We investigated potential effects of ecology on the developmental processes underlying ovariole number evolution among Hawaiian Drosophila, a large adaptive radiation wherein the highest and lowest ovariole numbers of the family have evolved within 25 million years. Previous studies proposed that ovariole number correlated with oviposition substrate [2-4] but sampled largely one clade of these flies and were limited by a provisional phylogeny and the available comparative methods. We test this hypothesis by applying phylogenetic modeling to an expanded sampling of ovariole numbers and substrate types and show support for these predictions across all major groups of Hawaiian Drosophila, wherein ovariole number variation is best explained by adaptation to specific substrates. Furthermore, we show that oviposition substrate evolution is linked to changes in the allometric relationship between body size and ovariole number. Finally, we provide evidence that the major changes in ovarian cell number that regulate D. melanogaster ovariole number also regulate ovariole number in Hawaiian drosophilids. Thus, we provide evidence that this remarkable adaptive radiation is linked to evolutionary changes in a key reproductive trait regulated at least partly by variation in the same developmental parameters that operate in the model species D. melanogaster.

Hawaiian picture-winged Drosophila exhibit adaptive population divergence along a narrow climatic gradient on Hawaii Island.

Anthropogenic influences on global processes and climatic conditions are increasingly affecting ecosystems throughout the world.Hawaii Island's native ecosystems are well studied and local long-term climatic trends well documented, making these ecosystems ideal for evaluating how native taxa may respond to a warming environment.This study documents adaptive divergence of populations of a Hawaiian picture-winged Drosophila, D. sproati, that are separated by only 7 km and 365 m in elevation.Representative laboratory populations show divergent behavioral and physiological responses to an experimental low-intensity increase in ambient temperature during maturation. The significant interaction of source population by temperature treatment for behavioral and physiological measurements indicates differential adaptation to temperature for the two populations.Significant differences in gene expression among males were mostly explained by the source population, with eleven genes in males also showing a significant interaction of source population by temperature treatment.The combined behavior, physiology, and gene expression differences between populations illustrate the potential for local adaptation to occur over a fine spatial scale and exemplify nuanced response to climate change.

Mapping Genomic Scaffolds to Chromosomes Using Laser Capture Microdissection in Application to Hawaiian Picture-Winged Drosophila.

Next-generation sequencing technologies have led to a decreased cost and an increased throughput in genome sequencing. Yet, many genome assemblies based on short sequencing reads have been assembled only to the scaffold level due to the lack of sufficient chromosome mapping information. Traditional ways of mapping scaffolds to chromosomes require a large amount of laboratory work and time to generate genetic and/or physical maps. To address this problem, we conducted a rapid technique which uses laser capture microdissection and enables mapping scaffolds of de novo genome assemblies directly to chromosomes in Hawaiian picture-winged Drosophila. We isolated and sequenced intact chromosome arms from larvae of D. differens. By mapping the reads of each chromosome to the recently assembled scaffolds from 3 Hawaiian picture-winged Drosophila species, at least 67% of the scaffolds were successfully assigned to chromosome arms. Even though the scaffolds are not ordered within a chromosome, the fast-generated chromosome information allows for chromosome-related analyses after genome assembling. We utilize this new information to test the faster-X evolution effect for the first time in these Hawaiian picture-winged Drosophila species.

A Test for Gene Flow among Sympatric and Allopatric Hawaiian Picture-Winged Drosophila.

The Hawaiian Drosophila are one of the most species-rich endemic groups in Hawaii and a spectacular example of adaptive radiation. Drosophila silvestris and D. heteroneura are two closely related picture-winged Drosophila species that occur sympatrically on Hawaii Island and are known to hybridize in nature, yet exhibit highly divergent behavioral and morphological traits driven largely through sexual selection. Their closest-related allopatric species, D. planitibia from Maui, exhibits hybrid male sterility and reduced behavioral reproductive isolation when crossed experimentally with D. silvestris or D. heteroneura. A modified four-taxon test for gene flow was applied to recently obtained genomes of the three Hawaiian Drosophila species. The analysis indicates recent gene flow in sympatry, but also, although less extensive, between allopatric species. This study underscores the prevalence of gene flow, even in taxonomic groups considered classic examples of allopatric speciation on islands. The potential confounding effects of gene flow in phylogenetic and population genetics inference are discussed, as well as the implications for conservation.

Incipient ecological speciation between successional varieties of a dominant tree involves intrinsic postzygotic isolating barriers.

Whereas disruptive selection imposed by heterogeneous environments can lead to the evolution of extrinsic isolating barriers between diverging populations, the evolution of intrinsic postzygotic barriers through divergent selection is less certain. Long-lived species such as trees may be especially slow to evolve intrinsic isolating barriers. We examined postpollination reproductive isolating barriers below the species boundary, in an ephemeral hybrid zone between two successional varieties of the landscape-dominant Hawaiian tree, Metrosideros polymorpha, on volcanically active Hawai'i Island. These archipelago-wide sympatric varieties show the weakest neutral genetic divergence of any taxon pair on Hawai'i Island but significant morphological and ecological differentiation consistent with adaptation to new and old lava flows. Cross-fertility between varieties was high and included heterosis of F1 hybrids at the seed germination stage, consistent with a substantial genetic load apparent within varieties through low self-fertility and a lack of self-pollen discrimination. However, a partial, but significant, barrier was observed in the form of reduced female and male fertility of hybrids, especially backcross hybrids, consistent with the accumulation of genetic incompatibilities between varieties. These results suggest that partial intrinsic postzygotic barriers can arise through disruptive selection acting on large, hybridizing populations of a long-lived species.

Deciphering the Routes of invasion of Drosophila suzukii by Means of ABC Random Forest.

Deciphering invasion routes from molecular data is crucial to understanding biological invasions, including identifying bottlenecks in population size and admixture among distinct populations. Here, we unravel the invasion routes of the invasive pest Drosophila suzukii using a multi-locus microsatellite dataset (25 loci on 23 worldwide sampling locations). To do this, we use approximate Bayesian computation (ABC), which has improved the reconstruction of invasion routes, but can be computationally expensive. We use our study to illustrate the use of a new, more efficient, ABC method, ABC random forest (ABC-RF) and compare it to a standard ABC method (ABC-LDA). We find that Japan emerges as the most probable source of the earliest recorded invasion into Hawaii. Southeast China and Hawaii together are the most probable sources of populations in western North America, which then in turn served as sources for those in eastern North America. European populations are genetically more homogeneous than North American populations, and their most probable source is northeast China, with evidence of limited gene flow from the eastern US as well. All introduced populations passed through bottlenecks, and analyses reveal five distinct admixture events. These findings can inform hypotheses concerning how this species evolved between different and independent source and invasive populations. Methodological comparisons indicate that ABC-RF and ABC-LDA show concordant results if ABC-LDA is based on a large number of simulated datasets but that ABC-RF out-performs ABC-LDA when using a comparable and more manageable number of simulated datasets, especially when analyzing complex introduction scenarios.

Genomic Signatures of Speciation in Sympatric and Allopatric Hawaiian Picture-Winged Drosophila.

The Hawaiian archipelago provides a natural arena for understanding adaptive radiation and speciation. The Hawaiian Drosophila are one of the most diverse endemic groups in Hawaiì with up to 1,000 species. We sequenced and analyzed entire genomes of recently diverged species of Hawaiian picture-winged Drosophila, Drosophila silvestris and Drosophila heteroneura from Hawaiì Island, in comparison with Drosophila planitibia, their sister species from Maui, a neighboring island where a common ancestor of all three had likely occurred. Genome-wide single nucleotide polymorphism patterns suggest the more recent origin of D. silvestris and D. heteroneura, as well as a pervasive influence of positive selection on divergence of the three species, with the signatures of positive selection more prominent in sympatry than allopatry. Positively selected genes were significantly enriched for functional terms related to sensory detection and mating, suggesting that sexual selection played an important role in speciation of these species. In particular, sequence variation in Olfactory receptor and Gustatory receptor genes seems to play a major role in adaptive radiation in Hawaiian pictured-winged Drosophila.

Postzygotic barriers isolate sympatric species of Cyrtandra (Gesneriaceae) in Hawaiian montane forest understories.

PREMISE OF THE STUDY: Recent reviews of reproductive isolation (RI) in plants propose that boundaries between closely related species are maintained predominantly through prezygotic mechanisms. However, few experimental studies have explored how boundaries are maintained in long-lived species. Hawaiian Cyrtandra presents an intriguing challenge to our understanding of RI, as it comprises 60 shrub or small tree species that are almost exclusively restricted to wet forests, where sympatry of multiple species is common. METHODS: We assessed the relative strengths of pre- and postzygotic barriers among four species of Cyrtandra occurring at the extremes of the main Hawaiian Island's natural island-age gradient, Kaua'i (4.7 Myr) and Hawai'i Island (0.6 Myr), to contrast the strengths and stages of reproductive isolation among species at different stages of divergence. KEY RESULTS: A combination of F1 seed germination, F1 seedling survival, and F1 seedling growth isolated (61-91%) three of the species from sympatric relatives. In contrast, the fourth species was isolated (59%) from its sympatric relative through phenological differences alone. Significant postzygotic barriers in between-island crosses were also observed in one species. CONCLUSIONS: Results suggest that boundaries between sympatric Cyrtandra species in Hawaii are maintained predominantly through postzygotic barriers. Observations from between-island crosses indicate that postzygotic barriers can arise in allopatry, which may be important in the initial divergence of populations. Future studies of RI in Cyrtandra should include a broader range of species to determine if postzygotic isolating barriers are foremost in the maintenance of species boundaries in this large genus.

Rapid adaptive radiation and host plant conservation in the Hawaiian picture wing Drosophila (Diptera: Drosophilidae).

The Hawaiian picture wing Drosophila are a striking example of adaptive radiation in specialist saprophages on an island system. We use DNA sequences from five nuclear genes with a total of 4260 nucleotides to provide a comprehensive phylogeny and biogeographic analysis of 90 species in the Hawaiian Drosophila picture wing clade. The current analysis indicates that the evolution of the picture wing clade took place more recently than previously suggested. The relationships of several morphologically anomalous taxa are resolved with strong support. Biogeography and host plant analyses show two periods of rapid divergence occurred when Kauai and Oahu were the main high islands, indicating that a combination of complex topographical features of islands and development of novel host plant associations was key to the rapid diversification of these lineages. For the past 2 million years, host associations within lineages have been largely stable, and speciation has occurred primarily due to the establishment of populations on newer islands as they arose followed by divergence by isolation. The existence of several apparently relictual taxa suggests that extinction has also played a major role in assembly of the present Hawaiian Drosophila fauna.